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Hyper-fast interstellar travel... 


• Is there a way within the framework of physics such that one 
could cross any given cosmic distance in an arbitrarily short 
period of time, while never locally breaking the speed of light 
(11 th commandment)? 

SPACEWARPS 



Images courtesy NASA 



Inflation: Alcubierre Metric 1 



Warp Metric: 
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Alcubierre, M., "The warp drive: hyper-fast travel within general relativity," 
Class. Quant. Grav. 11, L73-L77 (1994). 
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Bubble Topology Optimization 

York Time magnitude decreases 





bubble" thickness decreases 


Energy density magnitude decreases 


"bubble" thickness decreases 
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Surface plots of York Time & T°°, <v>=10c, 10 meter diameter volume, variable warp "bubble" thickness 
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Bubble Oscillation Optimization 
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Oscillate the bubble intensity 



Exotic Mass Warp Requirements, 10m diameter, v_ QDarent = 10c 


THINNER BUBBLE/RING 
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Energy Density Topology Approximations 


Chung -Freese null geodesics 
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Jupiter l.BxlO 77 
Earth 6.0xl0 M 
Moon 7.4X10 22 
4Vesta 2.6xl0 2D 
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Bulk Velocity, dU/dt 

■ dU/dt =0c 

* dU/dt=0.9c 

x dU/dt=0.999c 

* dU/dt =0. 99999c 

* dU/dt =0.99 999 99c 

+ dU/dt =0.99 999 99 99c 

Nimitz carrier 9.1xl0 7 

ISS 4.0x10 s 
Dragon^ 

^VoyagerT 




Warp Field Interferometer 



• Warp Field Interferometer developed after 
putting metric into canonical form 1 : 

ds 2 = (y 2 f(r s ) 2 -1)< dt — — dx > ~dx 2 +dy 2 + dz 2 

. /W -! , 

• Generate microscopic warp bubble that 
perturbs optical index by 1 part in 
10,000,000 

• Induce relative phase shift between split 
beams that should be detectable. 



Numerical simulation Numerical simulation 
of fringe at detector of fringe at detector 

with device off with device on 
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of contour p 
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White, H., "A Discussion on space-time metric engineering/' Gen. Rel. Grav. 35, 2025-2033 (2003). 





Warp Bubble Detection Attempts Using 

Interferometry 


* Goal is to use interferometer to detect and measure 
effect of warp bubble on optical path length through 
the measurement of associated interference fringe 
shifts. 

• The EWL has attempted to mitigate the effect of 
vibrations & air currents 

- Using a vibration-isolated optical table 

- Using a vibration-isolated room 

- Using an optical hood 

- Using signal processing to increase signal to noise ratio 

- Collecting statistical data to increase signal to noise ratio 
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Interferometer 
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Interferometer and Test-article Setup 




f-OVV-Vj 

Power 


f Computer for 
image capture 
article power c 


2 Polarizers 

(jointly used for 
intensity 
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Primary/Auxiliary 


Isolator pressure 


Height control pressure 
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SINGLE PATCH 

INTERFEROMETER MODAL ANALYSIS, 
ANALYTIC SIGNAL VARIATION 
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Interferometer Modal Analysis: 

Synthetic Data 

• A set of 1000 synthetic interferogram fringe patterns were generated as JPEG image 
frames 

- Set contained repeating sequence of 10 frames with sinusoidal shifting in time with a 
shifting amplitude of 1/10 of a spatial radian*. 

- No noise added to the shifting (i.e. unwanted vibrations not included) 

- Gray scale intensity sampled at upper left pixel and a Discrete Fourier Transform (DFT) 
calculated for the sample train. 

- Very strong peak observed at 1/5 Nyquist frequency as expected, and of course, no noise 
spectrum. 


shifted up and down along 
this axis 




2 fringe patterns out of set of 1000, showing maximum shift amplitude 
(In this case, the fringe pattern on the right is shifted about 1 radian up 
and to the left compared to the fringe pattern on the left.) 
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Modal Analysis: Analytic Signal Variation 

Single Patch 



The result is an array of complex values, called the analytic signal. Now, set all DFT coefficients to 

Use the complex value at the center and calculate the "phase". zero except for this small 

Repeat for all of the captured frames and store the representative ^ - region around this peak 

phase values as a ID series. corresponding to the fringe 

lines, and inverse transform^ 
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Single Patch Modal Analysis, 

always off 




Magnitude Spectrum of Phase Series 
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Single Patch Modal Analysis, 

cycled on/off 




•Test article cycled on and off at .25 Hz, and associated 
fringe images frame sampled at 2.5 Hz for 3149 samples. 
•Any effect should show up at frequency 315. 

•Nothing around that frequency stands out above the noise 
for this run. 

•No evidence that the test article caused fringes to shift at 
this measurement precision. 

•Concern that single patch may not sample enough of 
image, so subsequent effort done to sample across image. 





Magnitude Spectrum of Phase Series 
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MULTIPLE SAMPLE POSITIONS 

INTERFEROMETER MODAL ANALYSIS, 
ANALYTIC SIGNAL VARIATION 
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Description of Work 


• Room and table were floating 

• Optical hood covering table 

• Test article cycled on every 4 seconds 

• Frames captured every 0.4 seconds 

• 3149 frames processed 

• Sampling occurred at 88 positions in image plane 

• At each pixel sampled (through 3149 frames), the 
sample train was Fourier transformed, and the 
magnitude spectrum reviewed (see next slides) 
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Magnitude 


Test Article Off 



Magnitude Spectrum of Phase Series 



Frequency 



Magnitude 
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Test Article Cycled every 4 
(2 sec @ 20 KV, 2 sec off) 



Magnitude Spectrum of Phase Series 




Close up (off) 
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Close up (cycled) 
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Air Dielectric in Beam Path 

Cap ring replaced with high-voltage device "Ozona" such that 
the charge collection sphere partially blocked the interferometer 
beam in the leg along the laser axis. 

High voltage source was cycled on and off in 2 second intervals. 
The screen capture was operating at a rate of 0.4 sec per frame. 
The optics table was 
floating but not the 
room. 

The optics table was 
not covered. 

HV source power 
strip was on the 
micro flat table next 
to the interferometer 
table. 



°*«0 


HV Power Supply Off 
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HV Power Supply Cycled @ 2 second 

Intervals 
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DAN NEHLICH, SOUTH DAKOTA STATE UNIVERSITY 

INTERFEROMETRY MEASUREMENTS 
OF A DC TOROIDAL CAPACITOR 
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SDSU Interferometer Setup 





OS-*, 


f Image Subtraction Approach 

Positing that aberrations may include fringe pattern shifting in 
addition to or instead of induced curvature in the patterns, image 
subtraction was used to determine possible correlations to the 
voltage differential. 

- This method treats each image as an array and subtracts that array 
from the preceding, resulting in an image representative of the shift 
that occurred. 

- Average intensity of each subtracted image is calculated and statistical 
comparison is performed. 

- Statistical analysis of mean intensity values of each subtracted image, 
showed the observed shifts tended to increase during 
charge/discharge cycles, and remain constant in the control runs. 

The data presented in these scatter plots represents 2200 unique 
data points for charge/discharge cycles at 19 kV. 

Each point plotted represents approximately 20 images taken at the 
same point of the cycle, yielding an average intensity for each point 
in time along the graph. 



Control, Charge, Discharge Scatter Plots 



control 


charge 


discharge 


d(f) 

dt 
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Figure a: Image Subtraction during Test Article Discharge Cycle from 13 IsVto D V 







Hough Transform Approach 


In order to establish a direction of fringe 
pattern shifting, edge detection 
algorithms were necessary to determine 
individual line locations. 


The Hough Transform, along with the 
associated houghpeaks and houghlines 
functions in Matlab, was utilized to this 
end. 


Hough Line Representation 


An algorithm was written to utilize the 
edge mapping function of the transform 
determine the distance of each line to 
the origin (normal to the fringe 
orientation) and then sort and weight 
the data to provide the relevant fringe 
locations for analysis. 


Hough Transform with Peaks 


Average Shifts of Interference Pattern 


Cycle 

Average Shift 
(pixels/frame) 

St. Dev. 

Average Shift -limited 
( pixel s/frame) 

St. Dev. 

Control 

0.050652543 

1652178 

0.006819 178 

0 . 3 19178 

Charge 

-O.0128S0718 

1.720075 

0.010439022 

0.341535 

Discharge 

0.005357819 

O.S 10996 

0.015817141 

0.188S45 



ENHANCED WARP FIELD PHYSICS EXPERIMENTS 

IMAGE AVERAGING, FABRY-PEROT 
INTERFEROMETER, & TIME OF FLIGHT 
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VN 



Fabry-Perot Warp Field Interferometer 




Example: Michelson-Morley 
Interferometer image for 
Sodium source 



Consists of two reflecting, 
highly parallel surfaces, called 
an Etalon 

The interference pattern is 
created within the Etalon 

Multiple reflections in the 
Etalon reinforce the areas 
where constructive and 
destructive interference 
occurs 

Allows for much higher- 
precision measurements of 
fringes (image averaging 
without software) 


Example: Fabry-Perot 
Interferometer image for 
Sodium source (note doublet) 







SA-210 Scanning Fabry-Perot 
Interferometer 


Time of Flight Schematic 



200MHz test device 



Detectors Thorlabs High-Speed Avalanche Detector 



Model 

APD210 

Rise Time 

0.5 ns 

Supply Voltage 

+12 to +15 V 

Current Consumption 

200 mA 

Max. Incident Power 

10 mW 

Spectral Range 

400-1000 nm 

Frequency Range 

1-1600 mHz 

Maximum Gain 

2.5xl0 5 V/W 



HeNe 


EO Modulator 


Newport Broadband Amplitude Modulator 

Model 4102NF 

Type Broadband Amplitude Modulator 


Operating Frequency DC-200 MHz 


Wavelength Range 
Material 
Maximum Vn 


500-900 nm 
MgO:LiNbO s 
195 V @ 633 nm 


Maximum Input Power 2 W/mm 2 @ 532 nm 


Aperture Diameter 

2 mm 

RF Bandwidth 

200 MHz 

RF Connector 

SMA 

Input Impedance 

10 pF 

Maximum RF Power 

10W 

Connector 

SMA 




Agilent Technologies Infiniium DS09254A 2.5 GHz Oscilloscope 

•2.5 GHz bandwidth across all 4 analog channels 
•20 GSa/s max. sample rate 

•Standard 20 Mpts memory per channel, upgradeable to 1 Gpts 


36 





Forward Plan 


DC Toroidal Capacitor Approach: 

- Work with larger sample sets to decrease effects of vibrational "noise" and 
develop test articles with longer regions of optical influence to increase the 
signal. 

- Utilize image averaging algorithms to increase fidelity of interference 
information. 

- Utilize the Fabry-Perot Interferometer to increase the sensitivity of the 
experimental apparatus to below 1/100 of a wavelength. 


Explore the d<p/dt dependency in future test devices 

- The idea of an optimized space warp needs negative vacuum energy, and large 
d<f)/dt - both of these conditions are present in the q-thruster technology also 
being explored in the lab. 

- Use the q-thruster physics models to guide design of RF frequency test devices 
to be evaluated in the warp field interferometer, the Fabry-Perot 
Interferometer, and the time of flight experiment. 




Casimir force 




Quantum 

Vacuum 

MHD 


Q-th r usters are a low-TRL form of electric propu ision that operates on the 
principle of poshing off of the quant oro vacu u m . 

A terrestrial analog to this Is to consldler how a submarine oses Its propeller to 
posh a col u m n of water In one direction* while the sub recoils in the other to 
conserve momentum - the submarine does not carry a "ta n k" of sea water to he 
osed as propellant* 


In oor case* we ose the tools of Magnetohydrodynamlcs (MHD) to show how the 
thruster poshes off of the quantum vacuum which can he thought of as a sea of 
virtual particles - principally electrons and positrons that pop Into and out of 
existence* and where fields are stronger* there are more virtual particles. 

The Idea of poshing off the quantum vacuum has been In the technical literature 
for a few decades* hot to date* the obstacle has been the magnitude of the 
predicted thrust which has been derived analytically to be very small* and 
therefore not likely to be useful for human spaceflight. 

©or recent theoretical model development and test data suggests that we can 
greatly increase the magnitude of the negative pressure of the quantum vacuum 
and generate a specific force such that technology based on this approach can be 
competitive for in-space propulsion (''D.lN/kWj* and possibly for terrestrial 
applications (~10N/kW). 

As an additional validation of the approach * the theory allows calculation of 
physics constants horn first principles: Gravitational constant, Planck constant* 
Bohr radius* dark energy fraction* electron mass* 






C JS = U^,,. 
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Q-thruster Physics Data 


2004 Test Article 


jjjx 


2005 Test Campaign 


2012 Test Article 



ffl 


The test unit was run at 2.13 MHz, yielding an AC electric field 
of “20kV/m, and an AC magnetic field of “27 Gauss. 


• The observed thrust was +/- 0.89 mN 


The test unit was run at 3.8 MHz yielding an AC electric field 
of “20kV/m, and an AC magnetic field of “48 Gauss. 


• The observed thrust was +4.91 tc 
4900 mN (500gf) load cell 


-1.96 mN as measured 


As can be seen to the right, the thrust signal is very clear 
when the unit is excited. 


SFE Test Article at JSC 



In 2013, Boeing/DARPA sent Eagleworks Lab an SFE 
test article for testing and evaluation 

Evaluation of the test article in and out of a Faraday 
Shield performed from Feb through June 2013. 

• There is a consistent transient thrust at device 
turn-on and turn-off that is consistent with Q- 
thruster physics 

• The magnitude of the thrust scaled 
approximately with the cube of the input 
voltage (20-llQuN). 

• The magnitude of the thrust is dependent on 
the AC content of the turn-on and turn-off pulse 

® Specific force of transient thrust was in the ~1- 
20 N/kW range. 

-20-110 uN Thrust Pulses 
Specific Force ~1-20N/kW 



~3 mN Thrust 
Specific Force ~0.3N/kW 


Microwave Thruster Device 


SPR Ltd. Has produced several Microwave test articles. Claim is they produce 
thrust with just electric power input. 

• Shawyer's theoretical model has been deeme 

pwfoeWDcEnC witr.*!* i Ki Thruster assessed against Q-thruster models and 
analysis suggests this may be a microwave 
C version of a quantum vacuum plasma thruster. 

ir.. . • Tapered shape creates virtual toroid of active volume Prototype 1 6mN @ 850W, 

ST*. that can realize net thrust in virtual plasma. 0.02N/kW 

estricted to 


Unreal Results? 





- Thrust magnitude increased over multiple test 
devices from IS to 17QmN 
If Q-Thruster theory accounts for measured 
force, then microwave test articles may 
have ability to reach >10N/kW 
Chinese university claims to have duplicated 
EM Drive tests , but no way for U.S. to 
evaluate credibility (so we have Ignored it) 

16-170 mN Thruster 
Specific Force 0.02-0.4N/kW 



What we suspect: 

• A variety of industry experiments, for which theory is lacking, may be Q-thrusters 
including Boeing, Lockheed-Martin, EM Drive, Cannae, etc. 

• Low measured thrust but specific power ranges from 0.3 to 10+ N/KW 
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Q-thruster Roadmap 



s Pacelift ( 2oo, 


In 

-space 


Space-lift 



ero 

vm 






Q-THRUSTER + 2MW NUCLEAR POWER IS MISSION ENABLING 

POSSIBLE MISSIONS TO MARS, THE 
OUTER SOLAR SYSTEM, AND BEYOND 
WITH Q-THRUSTERS 
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Mars 




90 day stay 



• 140 day mission with 90 day 

stay at Mars 
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0.4 N/kW 
194 days 


90t spacecraft 

- 50t cargo, 20t power, 
20t propulsion 


• 2MW power 


JUPITER 


piter 
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i 


0.4 N/kW 
263 days 


SATURN 


SATURN 


4 N/kW 
86 days 


• 

90t spacecraft 

J^^ARTH 

• 

- 50t cargo, 20t power, 
20t propulsion 

2MW power 
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0.4 N/kW 
399 days 


JJRANUS 



4 N/kW 
129 days 


90t spacecraft 

- 50t cargo, 20t power, 
20t propulsion 


• 2MW power 
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Neptune 




APJH • 90t spacecraft 

- 50t cargo, 20t power, 
20t propulsion 

• 2MW power 


NEPTUNE 


NEPTUNE 


0.4 N/kW 
492 days 


4 N/kW 
160 days 
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Pluto 



PLUTO 



0.4 N/kW 
518 days 


4 N/kW 
167 days 


• 90t spacecraft 

- 50t cargo, 20t power, 
20t propulsion 


• 2MW power 
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0.4 N/kW 
5.6 years 


4 N/kW 
1.8 years 


90t spacecraft 

] - 50t cargo, 20t power, 
20t propulsion 

2MW power 



0.4 N/kW 
122.5 years 


CENTAURI 



4 N/kW 
29.9 years 


• 90t spacecraft 

- 50t cargo, 20t power, 
20t propulsion 

• 2MW power 
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Original Matthew Jeffries concept from mid 
1960% rendered by Mark Rademaker 



Rademaker\drex files 


Matthew Jeffries is the artist that created the familiar Star Trek enterprise 
look 
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Updated concept b 
by Mark R er^ 


Dr. White's theoretical findings, rendered 
|u/ork and inputs from Mike Okuda 
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Updated concept based on Dr. White's theoretical findings, rendered 
by Mark Rademaker with artwork and inputs from Mike Okuda 
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Reduced transit times for 


ISS Deployable impulsive trajectories 

or Class D Mission (Velocity Boost) 


Reduced low thrust 
transit time* 
(Reduced Alpha) 


1.5-50 Al 


0.1 AU 


TERRESTRIAL 
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• Interstellar Precursors 


Warp Field Interferometer 
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(Altered optics properties 
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Rapid Interstellar . 
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mission destinations 
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Principles of Q-thruster Operation 


Local mass concentrations, say in the form of a 
conventional capacitor with a ceramic dielectric, 
affect vacuum fluctuation density according to 
equation 1 


P 


v local 


= A 


P 


m local 


A 


= JP, 


m local P v 


(i) 


Just as relativistic acceleration (Unruh radiation) can 
change the apparent relative density of the vacuum, 
so too can higher order derivatives according to 
equation 2. 


Sp = 


8p = 


1 




4 ttG 

1 


1 ( da^\ 1 d 2 a 
+ 


V 


a 


V 


dt J 


ArTiG 


^ 1 ( d(f)^ 


\dt j 


V 





a = 


-V</> 


( 2 ) 


The tools of MagnetoHydroDynamics (MHD) can be 
used to model this modified vacuum fluctuation 
density analogous to how conventional forms of 
electric propulsion model propellant behavior. 





Gravitational Coupling Constant 


9 Consider the following thought experiment: what would an inertial observer in deep space find if 
the dark energy density were to be integrated over the light horizon of the observable universe, 
~13.7 billion light years? 

9 Starting with the Friedman Equation (and after some manipulation), the following equation can be 
derived that formally captures the results of this thought experiment: 


Vacuum Energy 
Density 




* Using 9.9xl0" 27 kg/m 3 [2] with t H of 13.7 billion years yields a predicted 
value for the gravitational constant of 6.45x10 11 m 3 /kg*s 2 

* A possible physical meaning to this rearranged equation solved for G is 
that gravitation is an emergent phenomenon rather than a fundamental 
force. 

* To be specific, the claim could be made that the gravitational coupling 
constant may be a long wavelength consequence (A =ct H ) of dark energy. 


Light horizon 
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Bohr Radius 


* The vacuum perturbation equation just derived can be used to evaluate the state of the 
quantum vacuum in close proximity of the proton at the center of the Hydrogen atom. 

* The first step is to calculate a quasi-classical density for the hydrogen nucleus. The radius of 
the hydrogen atom nucleus is given as /? 0 =1.2x10‘ 15 m {R=R 0 A' /3 where R 0 = 1.2x10‘ 15 m and A is 
the atomic number - these are experimentally determined by electron scattering). 

* The radius can be used with the mass of a proton to calculate a quasi-classical density of the 
hydrogen nucleus: 


• Using p^ 2 / 3 * 9.9x1 0 -27 kg/m 3 , along with this quasi-classical density p m , the perturbed negative 
pressure state of the quantum vacuum around the hydrogen nucleus is calculated to be: 


* The question can be asked how much volume of this perturbed state of the quantum vacuum 
is needed to have the equivalent energy value as the ground state of Hydrogen (13.6eV or 
2.1 8x1 O' 18 Nm) 


* The calculated radius is r= 5.29x1 0 11 m, which is an exact match to the given value for the 
Bohr Radius, a 0 = 5.29x1 O' 11 m. 






E 


r = 



Dr. Harold “Sonny” White 


Electron Mass 



Frank Wilczek, Nobel laureate: “We have achieved a beautiful and 
profound understanding of the origin of most of the mass of ordinary 
matter, but not of all of it. The value of the electron mass, in particular, 
remains deeply mysterious...” 

Consider the energy state of the perturbed quantum vacuum field around 
the proton, and set this equal to the kinetic energy of the orbiting electron 
at the ground state. 

4 3 2 1 2 

Pv local C =~ m e V 


We know the speed of the orbiting electron: 

v = ac = c / 1 37 

We can solve for the electron mass, and using the predicted value for 
Pv local °f 3.9x1 0 -5 kg/m 3 , we get a predicted electron mass of 9.1x10 -31 kg. 


8 


m e = 


m 0 Pv local C 


f \ 2 

C ' 


V 


137 


) 
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Magnetic Pressure 


* The first step now is to calculate the magnetic pressure around the Hydrogen 
nucleus. 


• The quasi-classical plasma pressure of the perturbed quantum vacuum state 
around the Hydrogen nucleus can be calculated by converting the electron velocity 
to temperature using 1/2 m e v 2 =kT, and making the assumption that the virtual 
electron-positron plasma has the same effective temperature as the orbiting 
electron. 

* When the plasma pressure calculation makes use of a 2/3 factor, analogous to the 
predicted dark energy fraction of 2/3 picked up during integration to calculate the 
Gravitational constant, the values are nearly identical: 


The magnetic field as perceived by 
the electron is given by the 
following relationship. The speed of 
the orbiting electron is ac. 



ju Q qv 


The magnetic pressure is a 


simple calculation: 

R 
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